Search for: All records

Global warming is expected to cause significant changes in the pattern of precipitation minus evaporation (P−E), which represents the net flux of water from the atmosphere to the surface or, equivalently, the convergence of moisture transport within the atmosphere. In most global climate model simulations, the pattern ofP−Echange resembles an amplification of the historical pattern—a tendency known as “wet gets wetter, dry gets drier.” However, models also predict significant departures from this approximation that are not well understood. Here, we introduce a new method of decomposing the pattern ofP−Echange into contributions from various dynamic and thermodynamic mechanisms and use it to investigate the response ofP−Eto global warming within the CESM1 Large Ensemble. In contrast to previous decompositions ofP−Echange, ours incorporates changes not only in the monthly means of atmospheric winds and moisture, but also in their temporal variability, allowing us to isolate the hydrologic impacts of changes in the mean circulation, transient eddies, relative humidity, and the spatial and temporal distributions of temperature. In general, we find that changes in the mean circulation primarily control theP−Eresponse in the tropics, while temperature changes dominate at higher latitudes. Although the relative importance of specific mechanisms varies by region, at the global scale departures from the wet-gets-wetter approximation over land are primarily due to changes in the temperature lapse rate, while changes in the mean circulation, relative humidity, and horizontal temperature gradients play a secondary role.

 

Abstract The response of zonal-mean precipitation minus evaporation ( P − E ) to global warming is investigated using a moist energy balance model (MEBM) with a simple Hadley cell parameterization. The MEBM accurately emulates zonal-mean P − E change simulated by a suite of global climate models (GCMs) under greenhouse gas forcing. The MEBM also accounts for most of the intermodel differences in GCM P − E change and better emulates GCM P − E change when compared to the “wet-gets-wetter, dry-gets-drier” thermodynamic mechanism. The intermodel spread in P − E change is attributed to intermodel differences in radiative feedbacks, which account for 60%–70% of the intermodel variance, with smaller contributions from radiative forcing and ocean heat uptake. Isolating the intermodel spread of feedbacks to specific regions shows that tropical feedbacks are the primary source of intermodel spread in zonal-mean P − E change. The ability of the MEBM to emulate GCM P − E change is further investigated using idealized feedback patterns. A less negative and narrowly peaked feedback pattern near the equator results in more atmospheric heating, which strengthens the Hadley cell circulation in the deep tropics through an enhanced poleward heat flux. This pattern also increases gross moist stability, which weakens the subtropical Hadley cell circulation. These two processes in unison increase P − E in the deep tropics, decrease P − E in the subtropics, and narrow the intertropical convergence zone. Additionally, a feedback pattern that produces polar-amplified warming partially reduces the poleward moisture flux by weakening the meridional temperature gradient. It is shown that changes to the Hadley cell circulation and the poleward moisture flux are crucial for understanding the pattern of GCM P − E change under warming. Significance Statement Changes to the hydrological cycle over the twenty-first century are predicted to impact ecosystems and socioeconomic activities throughout the world. While it is broadly expected that dry regions will get drier and wet regions will get wetter, the magnitude and spatial structure of these changes remains uncertain. In this study, we use an idealized climate model, which assumes how energy is transported in the atmosphere, to understand the processes setting the pattern of precipitation and evaporation under global warming. We first use the idealized climate model to explain why comprehensive climate models predict different changes to precipitation and evaporation across a range of latitudes. We show this arises primarily from climate feedbacks, which act to amplify or dampen the amount of warming. Ocean heat uptake and radiative forcing play secondary roles but can account for a significant amount of the uncertainty in regions where ocean circulation influences the rate of warming. We further show that uncertainty in tropical feedbacks (mainly from clouds) affects changes to the hydrological cycle across a range of latitudes. We then show how the pattern of climate feedbacks affects how the patterns of precipitation and evaporation respond to climate change through a set of idealized experiments. These results show how the pattern of climate feedbacks impacts tropical hydrological changes by affecting the strength of the Hadley circulation and polar hydrological changes by affecting the transport of moisture to the high latitudes. 

High‐resolution regional climate model (RCM) simulations of global warming consistently predict larger percentage increases in precipitation in the lee of midlatitude mountain ranges than on their windward slopes, indicating a weakening of the orographic rain shadow. This redistribution of precipitation could have profound consequences for water resources and ecosystems, but its underlying mechanisms are unknown. Here we show that rain‐shadow weakening is just one manifestation of a more general decrease in the influence of orography on precipitation under global warming. We introduce a simple model of precipitation change based on this principle, and find that it agrees well with an ensemble of high‐resolution simulations performed over the western United States. We argue that diminished orographic influence can be explained by the unique vertical structure of orographically forced ascent, which tends to maximize in the lower atmosphere where condensation is thermodynamically less sensitive to warming.

 

The influence of climate feedbacks on regional hydrological changes under warming is poorly understood. Here, a moist energy balance model (MEBM) with a Hadley Cell parameterization is used to isolate the influence of climate feedbacks on changes in zonal‐mean precipitation‐minus‐evaporation (P − E) under greenhouse‐gas forcing. It is shown that cloud feedbacks act to narrow bands of tropicalP − Eand increaseP − Ein the deep tropics. The surface‐albedo feedback shifts the location of maximum tropicalP − Eand increasesP − Ein the polar regions. The intermodel spread in theP − Echanges associated with feedbacks arises mainly from cloud feedbacks, with the lapse‐rate and surface‐albedo feedbacks playing important roles in the polar regions. TheP − Echange associated with cloud feedback locking in the MEBM is similar to that of a climate model with inactive cloud feedbacks. This work highlights the unique role that climate feedbacks play in causing deviations from the “wet‐gets‐wetter, dry‐gets‐drier” paradigm.

 

Abstract The stable isotope ratios of oxygen and hydrogen in polar ice cores are known to record environmental change, and they have been widely used as a paleothermometer. Although it is known to be a simplification, the relationship is often explained by invoking a single condensation pathway with progressive distillation to the temperature at the location of the ice core. In reality, the physical factors are complicated, and recent studies have identified robust aspects of the hydrologic cycle’s response to climate change that could influence the isotope-temperature relationship. In this study, we introduce a new zonal-mean isotope model derived from radiative transfer theory, and incorporate it into a recently developed moist energy balance climate model (MEBM), thus providing an internally consistent representation of the tight physical coupling between temperature, hydrology, and isotope ratios in the zonal-mean climate. The isotope model reproduces the observed pattern of meteoric δ 18 O in the modern climate, and allows us to evaluate the relative importance of different processes for the temporal correlation between δ 18 O and temperature at high latitudes. We find that the positive temporal correlation in polar ice cores is predominantly a result of suppressed high-latitude evaporation with cooling, rather than local temperature changes. The same mechanism also explains the difference in the strength of the isotope-temperature relationship between Greenland and Antarctica. 

This dataset contains monthly average output files from the iCAM6 simulations used in the manuscript "Enhancing understanding of the hydrological cycle via pairing of process-oriented and isotope ratio tracers," in review at the Journal of Advances in Modeling Earth Systems. A file corresponding to each of the tagged and isotopic variables used in this manuscript is included. Files are at 0.9° latitude x 1.25° longitude, and are in NetCDF format. Data from two simulations are included: 1) a simulation where the atmospheric model was "nudged" to ERA5 wind and surface pressure fields, by adding an additional tendency (see section 3.1 of associated manuscript), and 2) a simulation where the atmospheric state was allowed to freely evolve, using only boundary conditions imposed at the surface and top of atmosphere. Specific information about each of the variables provided is located in the "usage notes" section below. Associated article abstract: The hydrologic cycle couples the Earth's energy and carbon budgets through evaporation, moisture transport, and precipitation. Despite a wealth of observations and models, fundamental limitations remain in our capacity to deduce even the most basic properties of the hydrological cycle, including the spatial pattern of the residence time (RT) of water in the atmosphere and the mean distance traveled from evaporation sources to precipitation sinks. Meanwhile, geochemical tracers such as stable water isotope ratios provide a tool to probe hydrological processes, yet their interpretation remains equivocal despite several decades of use. As a result, there is a need for new mechanistic tools that link variations in water isotope ratios to underlying hydrological processes. Here we present a new suite of “process-oriented tags,” which we use to explicitly trace hydrological processes within the isotopically enabled Community Atmosphere Model, version 6 (iCAM6). Using these tags, we test the hypotheses that precipitation isotope ratios respond to parcel rainout, variations in atmospheric RT, and preserve information regarding meteorological conditions during evaporation. We present results for a historical simulation from 1980 to 2004, forced with winds from the ERA5 reanalysis. We find strong evidence that precipitation isotope ratios record information about atmospheric rainout and meteorological conditions during evaporation, but little evidence that precipitation isotope ratios vary with water vapor RT. These new tracer methods will enable more robust linkages between observations of isotope ratios in the modern hydrologic cycle or proxies of past terrestrial environments and the environmental processes underlying these observations.   Details about the simulation setup can be found in section 3 of the associated open-source manuscript, "Enhancing understanding of the hydrological cycle via pairing of process‐oriented and isotope ratio tracers." In brief, we conducted two simulations of the atmosphere from 1980-2004 using the isotope-enabled version of the Community Atmosphere Model 6 (iCAM6) at 0.9x1.25° horizontal resolution, and with 30 vertical hybrid layers spanning from the surface to ~3 hPa. In the first simulation, wind and surface pressure fields were "nudged" toward the ERA5 reanalysis dataset by adding a nudging tendency, preventing the model from diverging from observed/reanalysis wind fields. In the second simulation, no additional nudging tendency was included, and the model was allowed to evolve 'freely' with only boundary conditions provided at the top (e.g., incoming solar radiation) and bottom (e.g., observed sea surface temperatures) of the model. In addition to the isotopic variables, our simulation included a suite of 'process-oriented tracers,' which we describe in section 2 of the manuscript. These variables are meant to track a property of water associated with evaporation, condensation, or atmospheric transport. Metadata are provided about each of the files below; moreover, since the attached files are NetCDF data - this information is also provided with the data files. NetCDF metadata can be accessed using standard tools (e.g., ncdump). Each file has 4 variables: the tagged quantity, and the associated coordinate variables (time, latitude, longitude). The latter three are identical across all files, only the tagged quantity changes. Twelve files are provided for the nudged simulation, and an additional three are provided for the free simulations: Nudged simulation files iCAM6_nudged_1980-2004_mon_RHevap: Mass-weighted mean evaporation source property: RH (%) with respect to surface temperature. iCAM6_nudged_1980-2004_mon_Tevap: Mass-weighted mean evaporation source property: surface temperature in Kelvin iCAM6_nudged_1980-2004_mon_Tcond: Mass-weighted mean condensation property: temperature (K) iCAM6_nudged_1980-2004_mon_columnQ: Total (vertically integrated) precipitable water (kg/m2).  Not a tagged quantity, but necessary to calculate depletion times in section 4.3 (e.g., Fig. 11 and 12). iCAM6_nudged_1980-2004_mon_d18O: Precipitation d18O (‰ VSMOW) iCAM6_nudged_1980-2004_mon_d18Oevap_0: Mass-weighted mean evaporation source property - d18O of the evaporative flux (e.g., the 'initial' isotope ratio prior to condensation), (‰ VSMOW) iCAM6_nudged_1980-2004_mon_dxs: Precipitation deuterium excess (‰ VSMOW) - note that precipitation d2H can be calculated from this file and the precipitation d18O as d2H = d-excess - 8*d18O. iCAM6_nudged_1980-2004_mon_dexevap_0: Mass-weighted mean evaporation source property - deuterium excess of the evaporative flux iCAM6_nudged_1980-2004_mon_lnf: Integrated property - ln(f) calculated from the constant-fractionation d18O tracer (see section 3.2). iCAM6_nudged_1980-2004_mon_precip: Total precipitation rate in m/s. Note there is an error in the metadata in this file - it is total precipitation, not just convective precipitation. iCAM6_nudged_1980-2004_mon_residencetime: Mean atmospheric water residence time (in days). iCAM6_nudged_1980-2004_mon_transportdistance: Mean atmospheric water transport distance (in km). Free simulation files iCAM6_free_1980-2004_mon_d18O: Precipitation d18O (‰ VSMOW) iCAM6_free_1980-2004_mon_dxs: Precipitation deuterium excess (‰ VSMOW) - note that precipitation d2H can be calculated from this file and the precipitation d18O as d2H = d-excess - 8*d18O. iCAM6_free_1980-2004_mon_precip: Total precipitation rate in m/s. Note there is an error in the metadata in this file - it is total precipitation, not just convective precipitation. 

Abstract. Changes in glacier length reflect the integrated response to local fluctuations in temperature and precipitation resulting from both external forcing (e.g., volcanic eruptions or anthropogenic CO2) and internal climate variability. In order to interpret the climate history reflected in the glacier moraine record, the influence of both sources of climate variability must therefore be considered. Here we study the last millennium of glacier-length variability across the globe using a simple dynamic glacier model, which we force with temperature and precipitation time series from a 13-member ensemble of simulations from a global climate model. The ensemble allows us to quantify the contributions to glacier-length variability from external forcing (given by the ensemble mean) and internal variability (given by the ensemble spread). Within this framework, we find that internal variability is the predominant source of length fluctuations for glaciers with a shorter response time (less than a few decades). However, for glaciers with longer response timescales (more than a few decades) external forcing has a greater influence than internal variability. We further find that external forcing also dominates when the response of glaciers from widely separated regions is averaged. Single-forcing simulations indicate that, for this climate model, most of the forced response over the last millennium, pre-anthropogenic warming, has been driven by global-scale temperature change associated with volcanic aerosols. 

The hydrologic cycle couples the Earth's energy and carbon budgets through evaporation, moisture transport, and precipitation. Despite a wealth of observations and models, fundamental limitations remain in our capacity to deduce even the most basic properties of the hydrological cycle, including the spatial pattern of the residence time (RT) of water in the atmosphere and the mean distance traveled from evaporation sources to precipitation sinks. Meanwhile, geochemical tracers such as stable water isotope ratios provide a tool to probe hydrological processes, yet their interpretation remains equivocal despite several decades of use. As a result, there is a need for new mechanistic tools that link variations in water isotope ratios to underlying hydrological processes. Here we present a new suite of “process‐oriented tags,” which we use to explicitly trace hydrological processes within the isotopically enabled Community Atmosphere Model, version 6 (iCAM6). Using these tags, we test the hypotheses that precipitation isotope ratios respond to parcel rainout, variations in atmospheric RT, and preserve information regarding meteorological conditions during evaporation. We present results for a historical simulation from 1980 to 2004, forced with winds from the ERA5 reanalysis. We find strong evidence that precipitation isotope ratios record information about atmospheric rainout and meteorological conditions during evaporation, but little evidence that precipitation isotope ratios vary with water vapor RT. These new tracer methods will enable more robust linkages between observations of isotope ratios in the modern hydrologic cycle or proxies of past terrestrial environments and the environmental processes underlying these observations.